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We study a molecular engine constituted by a gas of A ~ 10^ molecules enclosed between a massive 
piston and a thermostat. The force acting on the piston and the temperature of the thermostat 
are cyclically changed with a finite period t. In the adiabatic limit t —^ oo, even for finite size 
N, the average work and heats reproduce the thermodynamic values, recovering the Carnot result 
for the efficiency. The system exhibits a stall time t* where net work is zero: for r < t* it 
consumes work instead of producing it, acting as a refrigerator or as a heat sink. At r > r* the 
efhciency at maximum power is close to the Curzorn-Ahlborn limit. The fluctuations of work and 
heat display approximatively a Gaussian behavior. Based upon kinetic theory, we develop a three- 
variables Langevin model where the piston’s position and velocity are linearly coupled together 
with the internal energy of the gas. The model reproduces many of the system’s features, such as 
the inversion of the work’s sign, the efhciency at maximum power and the approximate shape of 
huctuations. A further simplihcation in the model allows to compute analytically the average work, 
explaining its non-trivial dependence on r. 

PACS numbers: 05.70.Ln,05.40.-a,05.20.-y 


I. INTRODUCTION 

The usual statistical mechanics treats macroscopic ob¬ 
jects containing an enormous number N of particles, at 
least 0(10^°); the classical thermodynamics refers to adi¬ 
abatic processes. In practice a transformation can be 
considered adiabatic if its typical times are much longer 
than the times involved in the dynamics of the underlying 
system. Basically in the standard statistical mechanics 
and thermodynamics two asymptotic limits are present: 
large N and very slow changes of parameters [3,11. The 
challenge we face nowadays is going beyond these limits, 
extending thermodynamics and statistical mechanics to 
new models and applications Q. 

In fact, on one hand, it is clear that real transforma¬ 
tions occur in finite time: this problem has been fre¬ 
quently discussed in the recent past, giving birth to the 
so-called finite time thermodynamics [3, Hi, which fo¬ 
cuses on the study of engines working at finite power, 
i.e. far from Carnot efficiency. On the other hand the 
recent technological progresses allows us to relax also the 
large N limit: now it is possible to manipulate even small 
systems (say few hundreds particles) with non adiabatic 
changes of the parameters Q . Therefore it is necessary to 
(re)consider in details some aspects which are not partic¬ 
ularly relevant for macroscopic bodies. As an important 
example we mention the progresses in the study of fluc¬ 
tuations and their relation with response functions [3| • 

For the ambitious project of establishing a suitable 
statistical mechanics (as well as thermodynamics) for¬ 
malism for small systems and non adiabatic processes, it 
is necessary to build a theoretical framework, with new 
paradigmatic models, able to give an efficient descrip¬ 
tion of the statistical features at the mesoscale. The 
prototype of such a description is the Langevin equa¬ 
tion, which is able to catch the behavior of a colloidal 


particle (an object between the microscopic realm and 
the macroscopic one). The original Langevin equation 
has been established with a clever combination of macro¬ 
scopic arguments (the Stokes law for the friction force) 
and the use of statistical properties (equipartition). Fol¬ 
lowing the Smoluchowski approach to the Brownian mo¬ 
tion, sometimes, it is possible to rationalize the building 
of a Langevin equation. For instance for a big intruder in 
a diluted gas, using the kinetic the(^, one can derive the 
precise shape of the friction force Q- A step forward in 
this direction is represented by stochastic thermodynam¬ 
ics, based on the idea that the thermodynamical concepts 
of work, heat and entro py c an consistently be extended 
to a single trajectory [^, [l^ . 

The aim of the present work is the study of the ther¬ 
modynamic properties (work, heat, efficiency) in small 
systems performing non adiabatic cycles. We focus on 
a model system composed of a gas of molecules, a ther¬ 
mostat and a moving piston |lll4l3l |. proposing a partic¬ 
ular transformation (Ericsson cycle). The model’s re¬ 
sults are compared with those obtained for a coarse¬ 
grained stochastic equation which describes the evolu¬ 
tion of slowly changing quantities, such as position and 
velocity of the piston and temperature of the gas. Such 
a Langevin equation is obtained by means of kinetic the¬ 
ory considerations, in the spirit of the Smoluchowski ap¬ 
proach. 

The article is organized as follows: in Sec. |n]the model 
and the engine protocol are presented; in Sec. IIIII the 
results of numerical MD simulations of the model are 
reported, focusing on the average work and heat, the ef¬ 
ficiency and the fluctuations of such quantities. In Sec. 
m the coarse-grained stochastic model is derived and 
compared with the original system: finally, in Sec. |V] 
a simpler model is derived from the previous one and, 
in such a context, some analytical results are obtained. 
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The details of the calculations are discussed in the Ap¬ 
pendices. 

II. THE MODEL 



the study of the system in the canonical ensemble by 
means of standard statistical mechanics [Hill reveals 
that {Xy^ = {N + l)To/F anda^ = {N + 1)T^IF'^. In 
addition, if we define the estimate of the instantaneous 
temperature of the gas, 

1 ^ 

i=l 

the ensemble average of this quantity reads (T)®'? = To 
and its variance is cr|. = 2 /N. 


Figure 1. Sketch of the piston model. A gas of particles is 
confined by a fixed wall (the thermal bath) and a moving wall 
(the piston) that is subject to a constant external force. 


We consider an ideal gas of N pointlike particles with 
mass m, position and momentum (v^ = p^/m), 
i = 1... N, enclosed in a box. One of the sides of the box 
is a piston of mass M and momentum P which moves in 
the X direction. Indipendently of the real dimensionality 
of the box, only the motion in the x direction is relevant, 
as the particles interact only with the piston (see Fig. 
[1]). An externally controlled force F{t) = —Fx acts on 
the piston. The full (one-dimensional) hamiltonian of the 
system reads 


JL „2 p2 
■W = V -p — 
^ 2m 2M 

i—1 


FX, 


( 1 ) 


with the additional constraints 0 < Xi < A, f = 1,..., TV, 
and A > 0. The collisions with the piston are assumed 
to be elastic 


V' = V + 


2m 


m + M 


{vi - V), 


/ 2M , ^ 


( 2 ) 


where u' and V are post-collisional velocities. The wall 
at X = 0 acts as a thermostat at the temperature To', a 
collision of a particle with the wall is equivalent to give 
a new velocity v' to the particle with probability density 


/ /\ / — ILL^ _ \ 

e 2^0 0(u), 

J- O 


(3) 


A. The engine protocol 

When the parameters F and To vary in time, mechan¬ 
ical work can be extracted from the system. In partic¬ 
ular if we identify the (thermodynamical) internal en¬ 
ergy of the system with the value of the hamiltonian, 
E{t) — 'H(x(t),t), and, in addition, we define the input 
power as 




x(t) 


(5) 


conservation of energy simply reads E{t) = Q{t) + W{t), 
where Q is the energy absorbed from the thermal wall. 
For an hamiltonian TL as in Eq. © one gets W = XF. 
Let us remark that this formula is different from the one 
obtained in standard thermodynamics, W = FX: this is 
due to the fact that we included the energy of the piston 
in the internal energy of the system (T^ . 


F 



where Q{v) is the Heaviside Theta. Let us note that the 
presence of the piston introduces an interaction among 
the gas particles: for this reason the dynamics of the 
system depends on the number of particles N. As re¬ 
ported in previous studies of systems including the inter¬ 
action between a piston and one or more gases [TsI - fT^ . 
a relevant parameter for the dynamics is Nm/M: in¬ 
teresting behaviors are typically observed for values of 
this parameters 0{1), as in our study. Hereafter we 
use arbitrary units in numerical simulations, and we put 
fcs = 1 for the Boltzmann factor. At fixed F and To, 


Figure 2. Graph of F and To as a function of time over a 
cycle period r. 

Here we adopt the following cyclical protocol to obtain 
a heat engine: the parameters vary periodically in time 
over a cycle of length t (see Fig. [5] and the inset of Fig. 
|3] for a visual explanation). If we set to = kr, with k 
integer, the cycle has the following form: 

I) At times t G [to,to + r/4]: isobaric compression 
{F{t)=FL andTo(t) =Tff-4(Tff-Ti)(t-to)/r), 
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II) At times t G [to + T/4,to + '7‘/2]: isothermal com¬ 
pression {F{t) = Fl+ A{Fh - FL)[t - {to + r/4)]/T 
and To{t) =Tl)] 

III) At times tG [^0+^/2, to+3r/4] : isobaric expansion 
{F{t) = Fh and To{t) =Tl + A{Th — Tl))^ — (to + 
T/2)]/r); 

IV) At times t G [to-l-3r/4, to-l-r]: isothermal expansion 
(F(t) = Fh - 4(A// - FL)[t- {to + 3 t / 4 )]/t and 
To{t) = Th). 

The cycle of length r is repeated a large number of times 
over a long trajectory. In view of a Langevin-like analy¬ 
sis (see below) this protocol (called second type Ericsson 
cycle) - which is thermostatted for the whole duration 
of the cycle - is simpler than the more classical Carnot 
cycle. A similar model has been studied in [2^ with the 
difference that the velocity of the piston is given and can¬ 
not fluctuate (infinite mass M limit). The average values 
of heat and work in each segment of the cycle can be de¬ 
termined in the adiabatic limit, by substituting at every 
time t the value of each variable with the equilibrium 
average: X{t) = and E{t) = 

(see Table m. In particular, let us note that on segments 


Segment 

{W) 

{Q) 

I) 

0 

I{N+1){Tl-Th) 

11) 

(A+l)rUn(^) 

-{W) 

III) 

0 

^N+1){Th-Tl) 

IV) 

-(A+l)r«ln(^) 

-{W) 


Table I. Table with the adiabatic values of Q and W in each 
segment of the Ericsson cycle. The average (•) is intended 
over many realization of the cycle. 


II) and IV) no work is performed and that the heats ex¬ 
changed have same magnitude but opposite sign. There¬ 
fore, in the adiabatic limit, the system does not exchange 
net heat with any of the intermediate reservoirs at tem¬ 
perature Tl < T* < Th- Therefore in the rest of the 
paper we assume that two isobaric transformations do 
not contribute to the net exchange of heat and work: 
this is true for r —> oo and seems reasonable, for reasons 
of symmetry, at large r, while (small) discrepancies at 
finite r are observed in the simulations. We will denote 
with Qi the heat exchanged with the cold reservoir Tl 
in sector II), and with Q 2 the heat exchanged with the 
thermostat at temperature Th in sector IV). If (Q 2 ) > 0 
and {W) < 0, efficiency can be defined as 


{W) 






(Q 2 )’ 


( 6 ) 


where W is the total work on a cycle, and (•) denotes 
the average over many realizations of the cycle. Let us 
remark that this quantity is different from the average 
over many cycle of the fluctuating efficiency f] = WjQ^- 



Figure 3. Average values per cycle of work W and heats Qi, 
Q 2 as a function of the cycle typical time r. Dashed horizontal 
lines represent the adiabatic value of such quantities. Inset: 
schematic of the cycle protocol in the space of parameters 
F,To. 


III. RESULTS OF MD SIMULATIONS 


In order to perform molecular dynamics simulation of 
the system with time-dependent parameters F and To, we 
introduce an interaction potential between the piston and 
the particles to reproduce the effect of elastic collisions. 
We choose a repulsive soft sphere potential with cut-off 
radius tq: 


V{r) 



+ 12— - 13 
’"0 


0(ro - r), 


(7) 


where 0(r) is the Heaviside Theta. The parameter ro is 
to be chosen as small as possible, compatibly with inte¬ 
gration time-step At, in order to simulate a contact inter¬ 
action. In our case ro = 0.2 and At = 0.0005. The values 
of the other parameters, if not explicitly mentioned, are 
m = l,M = 100, N = 500, T^ = 11, Th = 13, Fl = 180, 
Fh = 220. The integration scheme adopted is based on 
the standard Verlet algorithm. 


A. Average work and heats 

In Fig. [3] we report the average values (over 500 cycles) 
of W, Qi and Q 2 as a function of r. In the adiabatic 
limit, T 10^, we recover, for {W), (Qi) and (Q 2 ), the 
values computed assuming quasi-static transformations 
in thermodynamics. At finite values of the cycle’s du¬ 
ration r < 10^ quite a complex scenario emerges. The 
absolute value of {W) decreases upon reducing t, until 
it vanishes at a stall time t* ^ 150. For shorter cycles, 
the engine consumes work instead of producing it (the 
regime at r > r* is marked, on Fig. [31 as “E”=engine). 
At smaller r, the analysis of the heats reveals the exis¬ 
tence of three regimes, marked on the Figure as “D”, “R” 
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and again “D”. In the “R” regime the system acts as a 
refrigerator, i.e. consumes work to push heat from to 
Th- In the “D” phases, the heat flow is the standard one 
(from Th to Tl), even if work is consumed: however, the 
rate of heat transfer {Q 2 )/t is higher than in the “E” 
phase, and therefore the machine acts as a more efficient 
heat sink, similar to dissipating fans. At a time r^es < t* 
we notice the presence of a maximum in {W): it is of the 
order of magnitude of the adiabatic limit, but with op¬ 
posite sign. At smaller r —^ 0 the consumed work goes to 
0. Let us note that the relevant timescales emerged from 
this analysis are in fair agreement with the characteristic 
relaxation times computed in a simple Langevin model 
of this system, see below. 


B. Power and efficiency 


Measures of the average developed power (P) = 
— {W)/t are reported in Fig. |4] (red curve) in the engine 
phase T > T* . Those measures are given as a function 
of the efficiency rj = —{W)/{Q 2 ), which is monotonically 
increasing with r. At a time around ^ 500, a maximum is 
observed in (P), whose associated efficiency is found to be 
slightly smaller than the Curzon-Ahlborn (CA) estimate 


r]cA = 1 — \/Tl/Th- We recall that the CA estimate 
is based upon an endo-reversible model of Carnot engine 
where the only entropy changes (even at finite times) oc¬ 
cur in the heat transfers. Recently a wider hypothesis 
for the CA formula has been investigated, which seems 
to be “symmetric dissipation”, i.e. equal entropy pro¬ 
duction rates during the two isothermals 22|. It is likely 
that our choice of values for Th and Tl = 0.85 Pf/ puts 
us close to that scenario. Nevertheless it is remarkable 
to recover a result usually derived through macroscopic 
theories, i.e. without fluctuations, in a small system such 
as our molecular model. 


C. Fluctuations 

In small systems, fluctuations are hardly negligi¬ 
ble [ 2 ^. In Fig. [5] (A and B, red curves), we display 
the behavior of fluctuations of work W integrated in a 
cycle for two different regimes, at t = 50 < r* and 
T = 500 > T*. Deviation from a Gaussian behavior are 
small, indicating that JV, even if finite, is large enough 
to expect the validity of the central limit theorem. In¬ 
terestingly the measure of the standard deviation (stdev) 
uw rescaled by the average value 5W = (Jw/\{W) \ (black 
curve in Fig. [5t)) shows that 5W <C 1 close to Tres and 
5W ^ I at the stall time r*. The relative stdev for the 
heat, 5 Q 2 = <^Q 2 /\{Q 2 )\ behaves much more regularly. It 
is also interesting to analyze the fluctuation of the “fluc¬ 
tuating efficiency”, i.e. fj = — WIQ 2 measured in a single 
cycle, see Fig. [Sp (restricted to positive values), which 
displays a long tail for values larger than the average [ 2 ^ . 
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Figure 4. Average power versus efficiency for MD (red curve) 
and for the reduced “3V” model, Eq. dill) (blue curve). The 
Curzon-Ahlborn estimate and the Carnot efficiency are also 




Figure 5. Study of fluctuations. A) and B): Pdf of the work 
in a cycle for two different values of r (“E” and “D” regimes), 
from the MD and from the reduced “3V” model Eq. (fTTI) . C) 
Pdf of the fluctuating efficiency in a cycle at r = 500. D) 
rescaled stdev (see text) of P{W) and of P(Q 2 ) as a function 
of r. The statistics in this Figure is obtained from 2000 cycles. 


IV. COARSE-GRAINED DESCRIPTION 


In order to make contact with stochastic thermody¬ 
namics [ 2 ^, which is a useful framework for small sys¬ 
tems, we need a coarse-grained description with few rel¬ 
evant (slowly-changing) variables. The contribution of 
the fast degrees of freedom is in the noise. Reasonable 
candidates are: the position of the piston X, its velocity 
V and the estimate of the instantaneous “temperature” 
of the gas T(t) = rn/N i^)■ The time evolution 

of these observables can be determined by computing the 
average rate of collision occurring between the particles 
of the gas and the walls of the container. Here, at any 
t, we assume the gas to be homogeneously distributed in 
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the interval ( 0 , X(t)) and each particle to have a velocity 
V, given by a Maxwell-Boltzmann distribution Pt(v) at 
the temperature T{t). In addition we use the fact that 
the collisions between the gas particles and the piston 
are elastic and that a particle that collides with the ther¬ 
mal wall gets a new velocity v' distributed according to 
a Maxwellian distribution po{v'). Taking into account 
the contributions of the external force and the collision, 
we have that the average derivative of the velocity of the 
piston iy) = limAt^o(^^)/^^ is 




F N 


&{v-V)pt{v). (8) 

(M + m) 


On the other hand, {T{t)) is the sum of two terms 
coming from the collisions with the piston 


(t) 


coll 


m 

X 


dv {v''^— v‘^)\v — V\Q{v — V)pt{v), (9) 


where v' is the velocity after an elastic collision, and the 
interaction of the gas with the thermostat 


(T) 


ther 


m 

X 


dv dv'{v'^ — v^) jt’l 0 (—v) pt(v)po(v'). 


( 10 ) 

In order to reduce the dynamics to a linear Langevin 
equation we assume the fluctuations of X, V and T to 
be small (such assumptions are reasonable if iV ^ 1 and 
M m) and expand Eqs. ([ 8 ]), (|9]) and (ITOl) up to the 
first order around the equilibrium values Xeq = NTq / F, 
Veq = 0 and Teg = To- The linearity of the equation is, 
on one hand, inspired by the gaussianity of pdfs, and, 
on the other, it is a useful assumption that allows simple 
computations. The stochastic part is obtained by adding 
the gaussian noise terms with amplitudes determined by 
imposing that the variances of the variables coincide with 
those computed within the canonical ensemble This 
yields 


V = -k(x - Xeq) --fVF yr - Teq) + 

2MTn 


2-fTo 

M 


Cl) 


T =-- 


N 


Ff,V-a{T-Teq) + ^J^^2, (II) 


where and ^2 are independent white noises (^i) = 0, 
mMt')) = kAt-t'), k{t) = F{t)y[MNTe{t)i 
7 (t) = 2F{t)y2m /\M^nTo{t)i p{t) = F(t) l[MTo{t)\ 
and a{t) = F{t)\j2/[mN’^iiTo(t)]. A numerical study of 
the “3 variables” (3V) model in Eq. (fTTl) reveals a fair 
agreement with our main observations. In Eig.[ 6 ]the aver¬ 
age values per cycle of work and heats are compared with 
those obtained in the original MD: we define W and Q 
as in MD, with E{t) = NT{t)/2 +MV{A/2 +F{t)X{t). 
The maximum and the inversion of the average work are 
fully reproduced, but with significant shifts of the valnes 
of r where they occur. Indeed, a more detailed analysis 
(not reported here) has identified the relevance of two 




Figure 6. Average work and heats per cycle. Comparison 
between the MD, the 3-variables model, Eq. (HU, and the 
analytical solution in Eq. m rescaled by their asymptotic 
values. 


additional variables: taking into account the position 
and the velocity of the center of mass of the gas, it is 
possible to achieve a better agreement with the MD. 
Unfortunately, the parameters for such a “5 variables” 
model can only be obtained by fitting the MD data. 
Notwithstanding its degree of approximation, the 3V 
model gives a fair account of efficiency at maximum 
power, see Eig. SI as well as of fluctuations, see Eig. 
and B, which have a very similar Gaussian shape and 
width. The overall shape of the efficiency fluctuations’ 
pdf (Fig. [5p) is also reprodnced. The eigenvalues of the 
dynamic’s matrix in Equations m give also access to 
typical timescales. For instance, for F = 200, T = 12, 
M = 100, m = 1 and N = 500, the eigenvalues read 
Ai ~ —0.02 and A 2,3 ~ —0.50 ± i 0.10 leading to three 
characteristic timescales to compare with the total dura¬ 
tion of the cylce: ti = 4/|Ai| ~ 170, T 2 = 4 /| 5 i(A 2 , 3 )| ~ 8 
and T 3 = 4 /| 5 (A 2 , 3 )| ~ 40. Remarkably, the order 
of magnitnde of the relevant timescales in the MD 
system is correctly reproduced by the eigenvalues of the 
equilibrium dynamic’s matrix (see Fig. |3] and Fig. [5]). 
A detailed stndy of the 3V model is out of our present 
scope, but certainly deserves future investigation. 


V. A SOLVABLE TOY MODEL 


In spite of its apparent simplicity, it is not easy to 
derive analytical results for the 3V model in a cycle of the 
external parameters. Here we show that the qualitative 
dependence of {W) on the total time of the cycle r can 
be obtained in a simplified version of Eq. CD, where we 
set the temperature T{t) to be equal to the temperature 
of the thermostat To(t) at every time t. In addition, 
we assume the parameters to vary in the following form 
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(w = 2-11 Jt) 

fit) = /o(l + ecos(a;i)), 

r(<) = ro(l + gesin(tLit)), (12) 

where f{t) = F{t)/M, e <C 1 and <7^0(1): we set 
/o = 2, To = 12, e = 0.1 and q = 0.8. In the adiabatic 
limit, this simplified cycle (an approximation of the Er¬ 
icsson protocol, see Fig. HD produces a work not very 
different from the one of the Ericsson cycle. Passing to 
average values {Y = (X)) we obtain the equation (see 
Appendix [BD: 


Y + k{t)Y + jit)Y = f{t). (13) 

The homogeneous solution associated to Eq. m goes 
to zero in the long time limit: therefore, since we are 
interested in the asymptotic stationary solution, we will 
focus only on the non-homogeneous solution. This will 
be done by expanding all the terms in Eq. (USD in powers 
of e. In particular, since Y{t) = Yo{t) + eYi{t) + 
by solving Eq. (IT^ for e —^ 0 one gets Yo{t) = and 

Yi + ujqYi -\- vYi = —/o cosLut + /ogsinwt, (14) 

where ujq = Fq /{MNTo) and v = 2fo\/2m/ {ttTo). The 
asymptotic solution is 

Yi (t) = A{uj) [cos(wt -I- 4>iuj)) — q sin(a;t -I- (('(w))], (15) 


A(w) = 


-/o 


(fiuj) = arctan 


(wq— w^) 

vui 


(16) 


The work performed over a cycle of the parameters of 
total time r can be now expressed in a simple way 

27r { t 

IT(r) = —M/o—e J dt (Fq + £^(0) sin ( 27r— 


= -M/ott a ( ^ 


sm ( 


27r\ , /27r 

— j ~ 9 ^ I — 


• (17) 


In Fig. [6] (black curve) it is seen that this result, when 
normalized to its adiabatic value, compares quite well, in 
spite of the many approximations introduced to obtain 
Eq. (HZD, with the average work performed by the MD 
system and the 3V model, recovering the change of sign at 
value not far from r* and a maximum at smaller values. 
Computing the heat transfers is a more difficult task, as 
it requires a solution of Eq. (IT^ at order O(e^). 


VI. CONCLUDING REMARKS 

Summarizing, we have introduced a new model for a 
heat engine where fluctuations (due to small X) and finite 


power (due to small r) are observed. In the past a great 
attention has been given to extremely simplified models, 
typically in the form of single molecules or colloids, or 
Markov chains inspired to biochemical reactions. Here 
we propose to move towards a higher level of complexity, 
and possibly realism, suggesting a new test-ground for 
statistical mechanics of small systems out-of-equilibrium. 
Our system reveals non-trivial features, such as several 
working regimes (engine, refrigerator, heat sink) tuned 
by simply controlling r. Notwithstanding its rich phe¬ 
nomenology, the model admits a coarse-grained descrip¬ 
tion in terms of three linearly coupled Langevin equa¬ 
tions. Further investigation of this reduced model are in 
order, in particular of heat, work and efficiency fluctua¬ 
tions [^ . 


ACKNOWLEDGMENTS 

We acknowledge useful discussions with T. Sano and 
D. Villamaina. Our work is supported by the “Granular- 
Chaos” project, funded by the Italian MIUR under the 
FIRB-IDEAS grant number RBID08Z9JE. 


Appendix A: Details on the derivation of the 
Langevin Equation 


In order to get a linear Langevin equation from ki¬ 
netic theory we start from the conditional equilibrium 
distribution with fixed values of the macroscopic vari¬ 
ables X, V, T, and then determine the average number of 
particles that, in the unit time At, collide with the piston 
or with the thermal wall. Using Eq. m and @ , one can 
determine post-collisional velocities and, accordingly, the 
average change of V and T, over a time At. In the fol¬ 
lowing, to simplify the notation, we denote with (•) the 
conditional average {■\X,V,T). 

The equation for X reads X = V. The total average 
force acting on the piston due to collisions is: 


(AU) 


lim 

A^-o At 


coll 


N 

X 


dv ■ 


2 m 


M 


{v — V)\V — v\pTiv)Q{V — v) = 


2Nm 


{M + m)X V 2t:T Jy 
N 

(m -f M)X 


dv (v - = 

v' 


(T -f mU") erfc 


2mrU2 

-e 2 t 


(Al) 


where erfc(a;) = 2/-^/^ dtexp(—t^). To obtain the 
total force, a term —F/M must be added. The elastic 
collisions with the piston also affect the kinetic energy of 
the gas, through the term 
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(T) = / dv {v'^ — v^)\v — V\Q{v — V)pt{v) = 

coll y(. I 


m / m 


X V 2ttT „ 
2M 


dv 




V + 2 


M 


m + M 


(y — u) ) — V 


{v - V)e = 


(m + M)‘^X IV TT 


2mT 


{2T - MV^) + y {MT + mMV^ - 2mT) erfc ( J ] . (A2) 


Finally, the average change of temperature in a time interval At, due to the collision with the thermal wall is simply 
given by the term 


(T) =17 dvdv'{v''^ -v‘^)\v\Q{-v)pTiv)po{v') = 

ther Ji. j 




XT„ V 27rT 


dv 


dv' {y'"^ - v"^) vv'e ’"are = 


2 Vt{To-T) 


X 


(A3) 


The equilibrium value of A, V and T for which (A) = 
0, (V) = 0 and T = (T) 


(T) 


+ {T) 

= 0 are 


coll 

ther 


Xeq = 

F ’ 


(A4) 

Veq = 

0 , 


(AS) 

Teq = 

To, 


(A 6 ) 


where terms 0{m/M) are neglected. We can obtain a 
linear equation by expanding the expressions above up 
to first order around the equilibrium values: this can be 
done only if fluctuations are small, i.e. when A ^ 1 and 
M ^ m. The sum of Eq. m and —F/M yelds 


{V) = -kiX-X,g)--fV + piT-T,g), (A7) 

with k = F'^/MNTo, 7 = 2 F2m/{M'^' kTo) and p = 
F/{MTo) ■ Similarly the sum of Eq. (IA2I) and Eq. (jA3l) 
yelds 

2MT 

= -^/rE-a(r-T,,), (A 8 ) 

with a = Fy^2/(mN’^nTo). The coefficients k, 7 , p and 
a vary in time according to the time evolution of F and 
Tq. In order to take into account the fluctuations of this 
variables one must add three independent gaussian terms 
= 'Cl and Ct = C 21 with an appropriate 
weight matrix by with = 0,1,2. In this particular 
ca se the matrix is d iagona l, bu = bi, with 60 = 0, 61 = 

^2'^ay and 62 = \/2'^a\. The final form of the linear 

Langevin equation thus reads 

A = E, 

y = _fc(X - Ae,) - 71^ + m(T - Te,) + 

rP _ 2 MTo frr rj. ^ ^ ^ ^ a\ 

T— -- a(T - Teg) + Y C 2 - (A9) 


Let us note that this equation, with fixed F and To, sat¬ 
isfies detailed balance (^ . 


Appendix B: Details on the analytic solution of the 
toy model 


Let us consider Eq. (US: 

Y + kit)Y + jit)Y = fit), (Bl) 


where 

, . 2F(t) I 2m 

“My nToit) ’ 

kit) = 

^ ’ MNToit) ’ 


(B2) 


and fit) = Fit)/M. In the Ericsson cycle F and T 
depend on time in a too complicated manner in order 
to perform analytic calculations. Therefore, in order to 
obtain an explicit result, in the following we will assume 


fit) = /o(l + ecos(wt)), 

Toit) = ro(l -I- gesin(a;t)), (B3) 


where uj = 27r/r, e <C 1 and q ^ 0(1). We will now 
sketch the derivation of the non-homogeneous solution of 
Eq. (IBII) as an asymptotic expansion in e <C 1. Let us 
expand in power of e all the coefficients appearing in Eq. 
(IBII) up to 0(e) 

Mp 

kit) Pi ° [1-|-(2coscLit — gsina;t)e] (B4) 

NTq 

lit) 

Yit)piYoit)+eYiit). (B5) 







































Plugging these expressions into Eq. (IBII) . for e = 0 
one gets 


y. + ^y. + 2/„^r„ = A. m 

leading to 

Yoit) = (B7) 

At the following order 0{e), Eq. (IB1|) gives 

Yi + -j^{Yi + {2cosojt- qsmojt)Yo) + 

JVJo 

2 / 0 ^^^ (^1 + (coswt - I sin Ujt^ Eq) = 

= fo cos Lut (B8) 


or, if we plug the value of Yq, 

Yi + WoEl + vYi = — fo cos cut + /ogsinwt, (B9) 

where Wq = yjMf^ /{NTo) and v = 2foyj2m/{ ttTo). A 
solution of this equation can be found in the form 

Yi (t) = A(a;)(cos(a;t + (/'(w)) + q sin(a;t + 0(w))), (BIO) 

where 



(Bll) 



(j>(oj) = arctan ( - - ) . 

(B12) 

\^o-^ J 
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